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We argue that the existence of Dark Matter (DM) is a possible consequence of GUT symmetry 
breaking. In GUTs like S'O(IO), discrete Z2 matter parity (— l) 3 ( s- - L ) survives despite of broken 
B — L, and group theory uniquely determines that the only possible ^2-odd matter multiplets belong 
to representation 16. We construct the minimal non-SUSY 50(10) model containing one scalar 16 
for DM and study its predictions below Mg. We find that EWSB occurs radiatively due to DM 
couplings to the SM Higgs boson. For thermal relic DM the mass range A/dm ~ 0(0.1 — 1) TeV is 
predicted by model perturbativity up to Ma- For Mom ~ O(l) TeV to explain the observed cosmic 
ray anomalies with DM decays, there exists a lower bound on the spin-independent direct detection 
cross section within the reach of planned experiments. 



I. INTRODUCTION 

The existence of Dark Matter (DM) of the Universe is now established without doubt [JJ. However, the 
fundamental physics behind it is unknown at present. In the most popular new physics scenario containing 
DM - supersymmetry - discrete i?-parity is imposed by hand to prevent phenomenological disasters such as 
fast proton decay [2]. Similarly, in dedicated DM extensions of the standard model (SM) with new singlet [3], 
doublet [1] or higher multiplct scalars [5], ad hoc Z 2 symmetry must be added to ensure the stability of DM. 
These phenomenological models cannot answer the two most fundamental questions related to DM: (i) why this 
particular multiplet or particle constitutes the DM of the Universe?; (ii) what is the origin of the imposed Z 2 
symmetry? Therefore the underlying physics principles related to the existence of DM remain obscured. 

In this work we argue that the existence of DM of the Universe can be a consequence of Grand Unification 
(GUT) . The GUT framework not only explains the origin of DM but also determines the type of the DM particle 
and constrains its properties. In this scenario the existence of DM, non-zero neutrino masses via seesaw [6] and 
baryon asymmetry of the Universe via leptogenesis [7] all point to the same GUT framework. 

We show that the Z 2 symmetry needed for DM stability could be a discrete remnant of GUT symmetry 
group, such as 5*0(10) [S] that we choose to work with in the following. When breaking 5*0(10) down to the 
SM gauge group SU(2)l X U(l)y, the 50(10) embedded U(l)x, where X is orthogonal to the SM hypercharge 
Y, leaves unbroken Z 2 [HI HU] 

Px = Pm = {-l? [B - L \ (1) 

which is the well known matter parity Pm- Due to its gauge origin Pm is a symmetry of any SM extension 
including non-SUSY ones. In the latter case group theory predicts uniquely, without any detailed model 
building, that the only possible Z2-odd multiplet under Eq. Q is the 16 of 50(10) [TTj . As inclusion of 
the fourth fermion generation I64 to the SM is not supported by experimental data, the non-SUSY 50(10) 
GUT predicts that the DM is a mixture of SU(2)l x U(1)y Pjvf-odd complex scalar singlet 5 and neutral 
component of doublet H 2 belonging to a new scalar 16 of 5O(10). Thus the DM of the Universe corresponds to 
the scalar analogues of the fermionic neutral matter fields, the right-handed neutrino Nr and the left-handed 
neutrino vl, respectively. Preserving Pm requires 5O(10) breaking by an order parameter carrying even charge 
of B — L [9j [10]. Therefore 5O(10) breaking also generates heavy Majorana masses which induce the seesaw 
mechanism as well as leptogenesis. 

To test the proposed DM scenario we study the scalar potential of a minimal 5O(10) GUT model containing 
one scalar 16 for the DM and one scalar 10 for the SM Higgs doublet. We derive [12] renormalization group 
equations (RGEs) for scalar mass parameters /if and interaction couplings Aj below the GUT breaking scale 
and study the vacuum stability and perturbativity conditions for those parameters. We find that the SM Higgs 
mass parameter n\ runs negative due to the presence of DM couplings with Higgs boson and triggers radiative 
electroweak symmetry breaking (EWSB) as in SUSY models [15] . Perturbativity up to Mq = 2 x 10 16 GeV 
restricts all scalar self couplings to be Ai < 1 at Mz predicting a restricted mass window 70 GeV % Mdm ^ 2 TeV 
for thermal relic DM mass. The operator mj ' {K^Mp)LLHiH 2 induces 3-body decays DM — > l~vW + which 
may explain the recently observed cosmic ray anomalies. For the DM mass preferred by this solution, Mdm ~ 
O(l) TeV, our framework predicts a lower bound on the spin independent direct cross section of DM with nuclei, 
which is within the reach of sensitivity of proposed experiments. 
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II. DM, LEPTOGENESIS AND SEESAW MECHANISM 

The 50(10) gauge group contains two orthogonal U(l) charges, which can be chosen to be the SM hypercharge 
Y remaining unbroken after 50(10) breaking at Mq, and broken X = 3(B — L) + 4T-sr, where T$r is the third 
component of SU(2)r £ 50(10) isospin. If 5O(10) is broken by fields with even X charge, the discrete subgroup 
Z 2 of U(l)x remains unbroken [9,- As 4T^r is always even, the surviving Px parity is nothing but the matter 
parity, Eq. ([IJ. Due to the 5O(10) breaking, B — L is broken at GUT scale generating large Majorana masses 
for right-handed neutrinos Nr. which suppress the light neutrino masses via the seesaw mechanism. The Nr. 
decays in early Universe induce the baryon asymmetry via leptogenesis. Thus, in our model, the existence of 
DM due to the matter parity Eq. , the existence of baryon asymmetry and the existence of seesaw suppressed 
masses of light neutrinos have the same GUT origin. To our knowledge, U(l)x, X = 5(B — L) — 2Y [T4] . 
has been used to forbid proton decay operators in GUTs, explicit examples of gauged U(1)b-l SUSY seesaw 
models generating i?-parity have been presented in [15) and low energy non-SUSY SM extension with extra 
U(l)' gauge symmetry generating Z 2 is presented in |16j . but the connection between non-supersymmetric DM 
and Px was first proposed in [TTj . 

The Px parity of 5O(10) matter multiplets is uniquely determined by group theory. Therefore the proposed 
5O(10) GUT scenario leaves no choice as to what are the DM particle multiplets. Under the group theoretic 
decomposition SU(5) x U(l) x the 16 representation of 5O(10) reads 16 = 1 16 (5) + 5 16 (-3) + 10 16 (1), where 
the X charges of the component multiplets are given in the brackets. While the X charges are different, all the 
fields in 16 of 5O(10) are odd under the conserved Z 2 parity Eq. pi. Interestingly, fields in 16 provide the only 
Px odd particles because all other fields coming from small 50(10) representations, 10, 45, 54, 120 and 126, 
are even under Px- Thus the SM fermions belonging to 16i, i = 1,2,3, are all P/^-odd while the SM Higgs 
boson doublet is P M -even because 10 = 5 10 (-2) + 5 10 (2). 

Although B — Lis broken in Nature by the heavy neutrino Majorana masses, discrete (— 1) 3 ( S ~ L ) is respected 
by the interactions of all matter fields. Therefore, without any model building, general GUT group theoretic 
argument implies that the non-supersymmetric DM must belong to 16 of 5O(10). Adding a new fermionic 16 
is equivalent to adding a new generation, which, due to mixing with lighter generations, cannot give DM. The 
only possibility is the new scalar 16 of 5O(10), which contains two DM candidates, the complex singlet 5 and 
the neutral component of the doublet H 2 ■ 



III. MINIMAL SO(10) GUT INDUCED DM MODEL 

The 5O(10) symmetric scalar potential of one 16 and one 10, 

V = (4 10 10 + Ai(10 10) 2 + nl 1616 + A 2 (1616) 2 
+ A 3 (10 10)(1616) + A 4 (16 10)(T610) 

+ X's [16 4 + h.c] + [16 10 16 + h.c] , (2) 

provides the minimal example of GUT DM model. Here we have taken all parameters to be real for simplicity. 
We assume 5O(10) to break at Mq down to SU(2)l x U(1)y X Pm in such a way that only one SM Higgs 
boson doublet H\ £ 10 and the DM candidate complex singlet 5 G 16 and the inert doublet H 2 £ 16 are light, 
with all other particle masses of order Mq. The 5O(10) symmetry breaking may occur in one or in several steps 
through intermediate symmetries such as SU(5) x U(l)x- We assume those steps to occur close to the GUT 
scale. Thus, between Mq and the EWSB scale M z , the DM is described by the Hx -> Hi , S -> -5, H 2 -> -H 2 
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invariant scalar potential 



V = ii{H\H t + X^Hl^f + \ilH\R 2 + X 2 (H J 2 H 2 f 

\2 



rt; 



rt; 



rt, 



+ Ms^S + [S 2 + (5t) 2 ] + Xs(S^S) 2 
+ >k [s 4 + (<?t) 4 ] + ^(St5) [S 2 + (St) 2 ] 
+ A S i(StS)(#Jfli) + A S2 (StS)(iftjj 2 ) 
+ ^i(ift^) [s 2 + (St) 2 ] + ^(fftff,) [s 2 



(St) 5 



(3) 



{H\H 2 f + {HlH,r 
S^H\H 2 +h 



XsiHjHJiHlHz) + A 4 ( J ff 1 t if 2 )(if 2 T ff 1 ) 

As 
2 



rt, 



SH[H 2 + h.c. 



together with the GUT scale boundary conditions 



rf(M G ) > 0, /i 2 (Af G ) = /4(M G ) > 0, 
A 2 (Af G ) = A S (M G ) = A S2 (M G ), A 3 (M G ) = A 5 i(M G ), 



(4) 



and 



Ml, ^<0(^)%?, 2 , (5) 

JW> J ^1.2,3,4- 

We require jU 2 (M G ) > in order not to break the SM gauge symmetry spontaneously at GUT scale. While the 
parameters in Eq. Q are allowed by SO (10), the ones in Eq. ^ can be generated only after SO(10) breaking 
by operators suppressed by n power of Planck scale Alp. If all parameters in Eq. (|5| vanished identically, Peccei- 
Quinn (PQ) symmetry would imply degenerate real and imaginary components of DM. However, direct search 
for inelastic DM requires the mass splitting to exceed 0(100) keV. Smallness of A5, as given by Eq. ([5]), allows 
one to interpret the annual modulation observed by DAMA experiment with inelastic scattering of DM in the 
Inert Doublet Model [17] . In our model there are more possibilities to obtain small mass splitting between real 
and imaginary components of DM candidates, cf. Eq. (|3|. In the following we assume the PQ symmetry to be 
broken softly by < < |/j 2 |. 



IV. VACUUM STABILITY CONSTRAINTS 



In the SM the requirements of vacuum stability and scalar potential perturbativity up to M G put the lower 
and upper bounds on the Higgs boson mass 127 GeV < Mh < 170 GeV, respectively (see [18] and references 
therein) . In our model the vacuum stability requires 

Ai > 0, 
A 2 > 0, 

Xs + X's>\^\, 
4A 1 (A S + A' S + A'^) > (A S1 + A^) 2 , 
4A 1 (A S + A' s -A'i) > (A S i-A^) 2 , 

Because there are more scalar couplings than in the SM, they can counteract the top quark Yukawa term in 
the beta function for Ai presented in the next Section and lower the vacuum stability bound on Mg below the 
LEP2 experimental bound of 114.4 GeV consistently with Eq. Therefore, in our model, the precision data 
indications for light SM Higgs boson does not contradict with vacuum stability constraints. 



A3 > — 2y A1A2, 

A 3 + A 4 - |A 5 | > -2VATA2", 

8(A S -A' S )A' S >A£ 2 , (6) 
4A 2 (A S + A' S + A'<0 > (A S2 + Ay 2 , 
4A 2 (A S + A' -A")>(A S2 -A') 2 . 
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V. RGE ANALYSES 



Because our low energy model is induced by GUT, it must stay perturbative up to the GUT scale. This 
requirement implies stringent constraints on the model parameters and consequently on the properties of DM. 
We have derived [T2J the full set of RGEs of the model Eq. The one-loop beta functions are given by 

(3 Xl = 2AX\ + 2A 2 + 2A 3 A 4 + A 2 + A 2 + X 2 S1 + 

+ j|(3.9 4 + 9' A + 2 3 V 2 ) - SA^g 2 + 5 ' 2 - 4yf) - 6yf, 
(3 X2 = 24A 2 . + 2A| + 2A 3 A 4 + A 2 + A 2 + A| 2 + A'| 2 

+ ^(35 4 + ff' 4 + 2 3 V 2 )-3A 2 (3 ff 2 + 5 ' 2 ), 
(3 Xs = 4(A X + A 2 )(3A 3 + A 4 ) + 4A 2 + 2A 2 + 2 A 2 + 2A S1 A S2 + 2A S1 A:< 



S2 

+ \{Zg* + g' 4 - 2g 2 g' 2 ) - 3A 3 (3g 2 + 5 ' 2 - 2y 2 ), 

/3a 4 = 4(Ai + A 2 )A 4 + 8A 3 A 4 + 4A 2 + 8A 2 

+ 3. 9 V 2 -3A 4 (3.g 2 + .g' 2 -22; 2 ), 
(3 Xb =4(A 1 + A 2 + 2A 3 + 3A 4 )A 5 

-3A 5 (3 5 2 + </ 2 -2y 2 ), 

p Xs = 20A 2 + 2A| 4 + Af x + 2A| 2 + A'| 2 + 36A'| + ^A s 2 , 

/3y s = A'# a + A'# 2 + 24A S A' S + 9 -Xf, 
% = 4A S iA' sl + 4A S2 A' S2 + 36(A S + \' s )\», 
pxsx = 4(3Ai + 2A S + Asi)A S i + 4A& + (4A 3 + 2A 4 )A S2 + 6A' S1 A S 

/3 Xs2 = 4(3A 2 + 2A S + A S2 )A S2 + 4A^ 2 2 + (4A 3 + 2A 4 )A S i + 6A S2 A S 

-l(3g 2 + g ,2 )Xs2, (7) 
Px> = (4A 3 + 2A 4 )A S2 + 4(3Ai + A s + 2A S1 + 3A S )A' S1 + 6A S iA s 



^(3.g 2 + 5 ' 2 -4 y 2 )A sl , 



tfg 2 + g' 2 )\' S2 . 



-§< 

/3a' S2 = (4A 3 + 2A 4 )A' S1 + 4(3A 2 + A s + 2A S2 + 3A S )A S2 + 6A S2 A S 

-I' 

h\ = 12 ^ Xl + 4 ^2^3 + 2m 2 A 4 + 2/i| Asi + 2/x' s 2 A' sl + - (/x| H + /x^) 

-^(3g 2 + .g' 2 -4y 2 ), 
^1 = 12^ 2 A 2 + 4 M 2 A 3 + 2 M 2 A 4 + 2/4 A S2 + 2// s 2 A S2 + i (/i| ff + /4 2 ff ) 

-^i(3 ff 2 + 3 ' 2 ), 

/V s = 8a4A s + 4 M 2 A 5 i + 4 M 2 A 52 + 6/is 2 A'^ + /4i? + Mshj 
/V s 2 = V#A S + 4/i 2 A' sl + 4/i 2 A' S2 + 12/^ 2 A' 5 + 6m|A s + 2fi S H»'sH, 
Push = 2^sif(A 3 + 2A 4 + Asi + A S2 ) + 2// SH (3A 5 + X' S1 + X' S2 ) 

3 '« 2 i „/2 o„,2\ 



^ SH (^ + g u -2yf i 



/V„„ = 2// SH (A 3 + 2A 4 + Asi + A S2 ) + 2fi SH (3X 5 + X' S1 + A' S2 ) 
'-u' SH (?>g 2 +g' 2 -2 y 2 ). 
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Figure 1: An example of Ai running from Mg to Mz- All \i not suppressed by SO(10) boundary conditions Eq. (J5J) are 
shown. 

We also include the one- loop /3-functions for g, g' ', (73 and yt, given by 



Based solely on the running due to the low-energy RGEs, we identify the unification scale 2 x 10 16 GeV by 
the RGE solution for g 2 = 53. The exact values of gauge couplings at Mq are given by g\ — yjhjig' — 0.58, 
92 = 93 — 0.53. Based solely on the running due to the low-energy RGEs, the precision of unification of all three 
gauge couplings in our model is better than in the SM because of the existence of an extra scalar doublet. We 
assume that an exact unification can be achieved due to the GUT thresholds corrections in full S'O(IO) theory 
which we cannot estimate because the details of GUT symmetry breaking are not known. Those corrections can 
have only small logarithmic influence on our numerical estimates of gi just below Mq and affect our numerical 
results for DM negligibly. In our numerical analysis we follow the strategy used in similar studies of parameter 
running in SUSY GUT theories. We fix all the measured model parameters and the SM Higgs boson mass at 
Mz and calculate the corresponding ix\ and \\. We run them up to the GUT scale where we randomly generate 
new physics parameters assuming the S*O(10) boundary conditions. We iterate running until the relative error 
between fixed and calculated n\ and Ai at Mz gets smaller than 1%. After that we calculate the DM abundance 
and direct detection cross section at Mz- 

We find that the new physics parameters are strongly constrained by the vacuum stability and perturbativity 
arguments. For example, assuming all A.; allowed by Eq. ^ to be equal at Mz, perturbativity of them up to 
Mq requires Xi(M z ) < 0.194, 0.187, 0.170 for M H = 120, 140, 160 GeV, respectively. We also impose the GUT 
boundary conditions Eq. Q . 

We present one consistent example of the evolution of scalar self couplings Ai and mass parameters /Uj between 
M z and M G in Fig. [I] and in Fig. [2] respectively. We assume M H = 140 GeV and take fj,' SH (M G ) = 1 GeV. 
The couplings Xi(Mz) must be small as not to reach the Landau pole below Mq. The SM gauge symmetry 
SU(2)l x U(1)y is not broken at Mq because all scalar mass parameters are positive ^(Mq) > 0. However, 
the SM Higgs mass parameter n\ exhibits stronger running than the DM mass parameters and triggers the 
radiative EWSB as in SUSY models [13]. In our case the EWSB is induced by DM couplings to the SM Higgs 
boson. (Previously, EWSB via a Coleman- Weinberg-likc mechanism has been considered in the Inert Doublet 
Model [T5].) 
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Figure 2: An example of /Ui,2,s running from Mg to Mz. Dashed line represents negative values of inducing EWSB. 

An interesting feature of the model, demonstrated in Fig. [2] is that the singlet DM mass parameter at low 
energies is always smaller than the doublet one, fig(Mz) < n\(Mz)- Thus, for small singlet-doublet mixing 
as is assumed in this example, the DM particle is predominantly scalar singlet S whose real and imaginary 
component mass degeneracy is lifted by small <§; |/i 2 |. 

VI. PREDICTIONS FOR DM MASS AND DIRECT DETECTION CROSS SECTION 

We assume that DM is a thermal relic and calculate its abundance and direct cross section with matter using 
the MicrOMEGAs package [20j. The DM interactions (pi were calculated using the FeynRules package [21] . We 
scan over the entire parameter space satisfying Eq. (RJ7 and calculate the RGE evolution of those parameters 
down to the EW scale. 

Fig. [3] presents a scattered plot of the spin-independent DM direct detection cross section per nucleon as 
a function of DM mass Mdm for the SM Higgs boson mass range from 115 GeV (red) to 170 GeV (violet). 
The whole parameter space allowed by theoretical constraints of vacuum stability and positive masses and 
experimental constraints from LEP2 and the WMAP 3ct result 0.094 < ^Inuh 2 < 0.129 [T] is shown. After fixing 
Ai and /i 2 from the assumed SM Higgs boson mass, we randomly generate the remaining scalar self couplings and 
mass parameters at Mq in the ranges < |Ai| < Air and < /if < (10 TeV) 2 . After RGE running the numerical 
ranges for non-zero parameters at M z are 0.117 < \ x < 0.239, 0.024 < A 2 < 0.227, -0.424 < A 3 < 0.247, 
-0.584 < A 4 < 0.599, 0.037 < A s < 0.177, 0.000 < \' s < 0.098, -0.212 < X S1 < 0.221, 0.031 < A S2 < 0.234, 
-14548 GeV 2 < [i\ < -7056 GeV 2 , 3147 GeV 2 < [4 < 1.72 x 10 7 GeV 2 , 2894 GeV 2 < m| < 3.84 x 10 6 GeV 2 , 
-11634 GeV < n' SH < 11504 GeV. 

The present [22] and future [23] experimental sensitivities for DM direct detection are shown in Fig. [3j As 
a result, we find that the DM mass is restricted to the window 70 GeV < Mdm ^ 2 TeV. The lower bound 
comes from non-observation of charged scalars at LEP2. The upper bound A/dm ^5 2 TeV comes from the the 
requirement of perturbativity of the model parameters up to Mq. Therefore, the DM mass scale Mdm ^ 0(0.1 — 
1) TeV is a prediction of our scalar DM GUT scenario. 

The direct DM interaction with nuclei occurs via the SM Higgs boson exchange. The dominant DM-Higgs 
effective coupling involved in this process is 



Acs v = -{V2scfj,' SH + 2s 2 (A 3 + A 4 )u + 2c 2 \ S iv), 



(9) 
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Figure 3: DM direct detection cross-section per nucleon vs. A/dm- Color shows SM Higgs masses from 115 GeV (red) 
to 170 GeV (violet). The points shown encompass the whole parameter space allowed by theoretical and experimental 
constraints. 

where s, c are the sine, cosine of the singlet-doublet mixing angle. If A/dm ^ 300 GeV, cancellation between 
different terms in Eq. (|9| is possible and the spin independent direct detection cross section can be accidentally 
small, cf. Fig. [3] However, for larger DM masses both Eq. ^ and thermal freeze-out cross section are dominated 
by large fi' SH term and one obtains a relation between the DM abundance and the direct detection cross sections 
with only mild dependence on Mh via RGEs. For Mdm = 1 TeV the WMAP result predicts a lower bound 
a/n > 2 • 10~ 43 (115 GeV /Mh) 4 cm 2 which is well within the reach of the planned experiments, cf. Fig. [3] 
Should CDMS experiment observe DM scattering on nuclei, light SM Higgs boson mass is indicated in agreement 
with precision electroweak data. 

VII. DM INDIRECT DETECTION 

The PAMELA [S3], ATIC 125], HESS [2E] and Fermi [27] anomalies of cosmic ray positron/electron fluxes 
can be explained with 0(1) TeV mass DM decays [28] via d = 6 operators [29], preferably to multi-particle final 
state [301 I3T] . Non-observation of photons associated with DM annihilation in the Galactic center [32] and in 
DM haloes in the Universe [33] as well as the suppression of hadronic annihilation modes [33] strongly favor 
DM decays over annihilations as a solution to the anomalies. 

In our scenario the decays of DM are most naturally explained via the seesaw like operator LLH1H2 which, 
in addition to the suppression by heavy Majorana neutrino scale Mjv, must be suppressed by the Z 2 breaking 
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effects by additional Mp. We obtain that below EWSB scale the dominant decay mode is given by 

m -LLHiH 2 -> 1CT 30 GeV"Vr W + H$, (10) 



M N M P 

where we have taken Ajv ~ 1, Mjv ~ 10 14 GeV and m ~ v ~ 100 GeV. In the decays of 1F + antiprotons are 
produced about in 10% of decays. Such a small fraction of antiprotons is till allowed by PAMELA data taking 
into account uncertainties in the cosmic ray propagation models |35j . Such a small effective Yukawa coupling of 



Eq. (10) can explain the long DM lifetime 10 s without conflicting with the present observational constraints. 



VIII. CONCLUSIONS 



We have argued that the existence of DM, the baryon asymmetry of the Universe, and small neutrino masses 
may all signal the same underlying GUT physics. Although B — L is broken in Nature by heavy neutrino 
Majorana masses, Z 2 parity (— 1) 3 ( S ~ L ) is respected by interactions of all matter fields. Hence, group theory 
predicts that in S'O(IO) GUTs the non-supersymmetric DM must be contained in the scalar representation 16. 

Based on S'O(IO) GUT, we have presented a minimal DM model, calculated the full set of its RGEs and 
studied its predictions. We find that the EWSB occurs radiatively due to SM Higgs boson couplings to the DM, 
analogously to SUSY models. The thermal relic DM mass is predicted to be A/dm 0(0.1 — 1) TeV by the 
requirement of perturbativity of model parameters up to the GUT scale. If M u m > 300 GeV as suggested by 
DM decay solution to the recently observed cosmic ray anomalies, the WMAP measurement of DM abundance 
predicts a lower bound on DM spin independent direct cross section with nuclei, which is within the reach of 
planned experiments for all values of the SM Higgs boson mass. 
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